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ABSTRACT 

We compare galaxy number counts in Advanced Camera for Surveys (ACS) fields containing 
moderate-redshift (0.2 < z < 1.0) strong gravitational lenses with those in two comparison 
samples: (1) the first square degree of the COSMOS survey, comprising 259 ACS fields and 
(2) 20 "pure parallel" fields randomly located on the sky. Through a Bayesian analysis we 
determine the expectation values (fio) and confidence levels of the underlying number counts 
for a range of apertures and magnitude bins. Our analysis has produced the following results: 
(i) We infer that our control samples are not consistent at the >lQ-a level, with the num- 
ber counts in the COSMOS sample being higher than in the pure parallel sample. This result 
matches those found in previous analyses of COSMOS data using different techniques, (ii) We 
find that small-size apertures, centered on strong lenses, are overdense around the l-a level 
compared with randomly placed apertures in the control samples, even compared to the COS- 
MOS sample. Correcting for the local clustering of elliptical galaxies, based on the average 
two-point correlation function, this over density reduces to the < l-cr level. Thus, the over- 
density of galaxies seen along a typical line of sight to a lens can be explained by the natural 
clustering of galaxies, rather than being due to lenses lying along otherwise biased lines of 
sight, (iii) Despite the considerable scatter in the lines of sight to individual lens systems, we 
find that quantities that are linearly dependent on the external convergence (e.g. Hq) should 
become unbiased if the few extra galaxies that cause the bias (i.e. AyUo ^ 2.0 galaxies with 
19 < m < 24 for aperture sizes < 45" radius) can be accounted for in the lens models. 

Key words: gravitational lensing - large-scale structure of Universe - distance scale 



1 INTRODUCTION 

Strong gravitational lenses, where multiple images of the back- 
ground object are formed, are powerful probes of the distribution 
of mass in the Universe. The properties of the lensed images are, in 
principle, sensitive only to the projected mass of the lensing object, 
with no requirements that the mass be lumin ous or baryonic (e.g., 
ISchneider. Kochanek. & Wambsganssll2006l) . Most of the lensing 
signal comes from the primary lensing object - typically a mas- 
sive early-type galaxy - and, if the lensing object is a member of a 
galaxy cluster or group, its immediate environment. 

However, the distribution of large-scale structure (LSS) 
along the line of sight to the lens system adds perturba- 
tions to the lensing properties. For example, simulations have 
shown that a non-negligible fraction of lenses can only be pro- 
duced by having multiple lens plan e s along the line of sigh t 
jWambsganss. Bode. & OstrikeJbOOSl : Iffilbert et alj|200l |2008|) . 
Furthermore, the differences in light travel times along rays form- 
ing the multiple lensed images in a given sys tem can be affected at 
the level of a fe w percent (e.g.. lSelia3l 19941) or up to ~10% (e.g., 
iBar-Kanal 19961) by the distribution of LSS along that particular line 
of sight. These effects should be random for random lines of sight. 



so it should be possible to reduce the LSS uncertainties and exploit 
the power of gravitational lenses as cosmological tools by averag- 
ing over many systems. If, however, lenses lie along biased lines of 
sight, this reduction will not occur and global parameters such as 
Ho determined from large lens samples will be biased. 

To date, most observational investigations of the effects of the 
environment on strong lensing have focused on the local neighbor- 
hood of the lens by searching for spectroscopic evidence of galaxy 



ing galaxv (e.g., jKundic et alj 1997allb 


; Fassnacht & LubinI 


20021; 


Momcheva et alJuOOa: iFassnacht et al 


l2006hl: lAuser et alj 


2007|, 


2008"). Some of these investigations (e.g.,'Fassnacht & Lubin' 


20021; 


Momcheva et al. 2006; Fassnacht et al. 2006b; Auger et al. 200% 



have also found mass concentrations along the line of sight that 
are at different redshifts than the lensing galaxy. However, due 
to the limitations imposed by requiring spectroscopic redshifts for 
their analyses, the spectroscopic surveys are necessarily incomplete 
samples of the line of sight. Most also are also biased because they 
preferentially target galaxies expected to be at the redshift of the 
primary lens. There have also been photometric studies of lens 
fields in order to evaluate lens environments. These photometric 
investigations are the closest in concept to the analysis in this pa- 
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per. However, t hey focused on either group finding via detection of 
red sequences jWilliams et al.ll20o3) . or on describing the imme- 
diate environment of the lensing galaxy by using galaxy c olors to 
strongly favor galaxies likely to be at the redshift of the lens jAugeJ 
l2008| - |Treu et al.l2009l) . None of these previous photometric studies 
evaluates the contributions by galaxies along the full line of sight. 

In this paper, we investigate the question of whether lens lines 
of sight are biased by comparing, through Bayesian and frequentist 
statistics, the number counts of galaxies in fields containing gravi- 
tational lenses with those obtained from two control samples. The 
control samples are chosen to provide reasonable approximations 
to typical lines of sight through the Universe. All image s were ob- 
tained with the Advanced Camera for Surveys (ACS: [pord et al.l 
Il998l l2002h aboard the Hubble Space Telescope (HST). The un- 
derlying idea is that lines of sight that are overdense in galaxies 
are also overdense in mass, if the underlying redshift distributions 
are roughly the same. Hence by simply counting galaxies, one can 
make conservative statements about the lines of sight towards lens 
galaxies that are not highly model-dependent. We discuss the sam- 
ple selection and image processing in §2, do a frequentist analysis 
of the samples in §3, develop and use a Bayesian framework for 
comparing the samples in §4, briefly describe the number counts 
for individual lens systems (as opposed to sample averages) in §5, 
and interpret the results in §6. 



Table 1. Lens Sample 









Lens System 


(sec) 


■2^1 ens 


JVAS B0218+357 


48720 


0.68 


CLASS B()445+128 


5228 


0.557 


CLASS B0850+054 


2296 


0.59 


CLASS B 1608+656 


2528 


0.630 


CLASS B2108+213 


2304 


0.365 


CFRS 03.1077 


2296 


0.938 


HE 0435-1223 


1445 


0.46 


HE 1113-0641 


1317 


0.75" 


J0743+1553 


2300 


0.19 


JOS 16+5003 


2440 




J 1004+ 1229 


2296 


0.95 


Q1131-1231 


1980 


0.295 


SDSS 0246-0825 


2288 


0.72a 


SDSS 0903+5028 


2444 


0.388 


SDSS 0924+0219 


1148 


0.39 


SDSS 1004+4112 


2025 


0.68 


SDSS 1138+0314 


2296 


0.45 


SDSS 1155+6346 


1788 


0.176 


SDSS 1226-0006 


2296 


0.52 


WFI 2033-4723 


2085 


0.66 



Redshifts marked with an " are photometric redshifts. All other redshifts 
are spectroscopic. 



2 DATA REDUCTION 

In this section we briefly describe the data reduction and catalogue 
extraction. 

2.1 Sample Definition and Data Processing 

The lens sample comprises 18 systems which were observed with 
ACS as part of the CfA-Arizona Space Telescope Lens Survey 
(CASTLES; GO-9744; PI Kochanek). This sample was defined by 
taking the full list of lenses observed with ACS as part of the CAS- 
TLES program (24 galaxies in all) and only including systems that 
(1) had total exposure times > 2000 s, (2) had no extremely bright 
stars in the field, and (3) had galactic latitudes of |6| > 10. Each 
system was observed through the F555W and F814W filters, but 
for comparison with the control data sets, we only consider the 
F814W data. The typical total exposure times through the F814W 
filter were ~2000-3000 sec. Details of the observations are given in 
Table I . The pipeline-processed data were obtained from the Multi- 
mission Archive at Space Telescope, and the ind ividual exposures 
were combined using the multidrizzle package jKoekemoer et al.l 
l2002h . We also included in the lens sample deep ACS images 
of B02 18+357 (GO-9450; PI N. Jackson) and B 1608+656 (GO- 
IOI58; PI Fassnacht), with total F814W exposure times of 48,720 
and 28,144 sec, respectively. For easier comparison with the rest of 
the lens sample we only used the data from the first four F8I4W ex- 
posures on the B 1608+656 field, with a combined integration time 
of 2528 sec. The lens galaxies are at moderate redshifts, with a 
mean and RMS of = 0.55 and — 0.22, respectivel y. This 
should be c ompared to the en vironmental investigations of lAugeil 
( |2008|) and lTreuetai] | |2009|) , which focused on a lower redshift 
sample (z < 0.3). 



2.2 Control fields 

The first control sample consists of data obtained by the Cosmic 
Evolution Survey team (COSMOS: [Scoville et al1l2007jlbh . The 
COSMOS data consist of a mosaic of approximately two square 
degrees, with all of the images obtained through the F814W fil- 
ter. They were obtained as part of a 510-orbit HST Treasury pro- 
posal in Cycles 12 and 13 (GO-9822; PI Scoville). Each field has 
a total exposure time of 2028 sec, comparable in depth to the lens 
fields. The data have b een fully reduced by the COSMOS team 
jKoekemoer et al.ll2007l) . so it was not necessary to run multidriz- 
zle. Instead the 257 processed science and weight images from Cy- 
cle 12 - comprising ~1 square degree - were obtained from the 
COSMOS ACS website hosted by the NASA/IPAC Infrared Sci- 
ence Archiv^ 

The second control sample consists of data obtained as part of 
a pure parallel program to search random fields for emission line 
galaxies (GO-9468; PI L. Yan). This program included 28 point- 
ings in the F814W filter, of which we used the 20 that success- 
fully passed all of the criteria required for our image processing 
(hereafter referred to as the "pure parallel fields"). These pure par- 
allel fields cover a total area of ~0.06 square degrees. Although the 
number of pointings is much smaller than obtained with the COS- 
MOS program, the pure parallel fields have the advantage of not 
being contiguous on the sky and, thus, provide an important check 
that the COSMOS data do not have some overall bias in the num- 
ber counts due to sample variance. These data were processed by 
the a modified version of the pipeline developed as part of the HST 
Archive Galaxy Gravitational Lens Survey (Marshall et al., in prep) 
which is designed to produce final images aligned sufficiently well 
to conduct weak lensing analyses. The modification to the pipeline 
for the pure parallel fields was simply to change the output pixel 
scale from the HAGGLES standard 0'.'03 pix"^ to 0'.'05 pix"\ in 
order to match the COSMOS pixel scale. 

^ http://irsa.ipac.caltech.edU/data/COSMOS/images/acs_vl.2/ 
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2.3 Object Detection and Flagging 

Object catalogs were obtained by running SExtractor 
teertin & ArnoutsI 1 19961) on each ACS image. For both the 
lens-field and control files, the weight maps produced by multidriz- 
zle were used to improve the object detection. The magnitudes 
used in this paper are Vega-based FHIAW magnitudes measured 
within the "AUTO" aperture computed by SExtractor. Hereafter, 
we will use m to designate these F814VF Vega magnitudes. We 
edited the initial catalogs to reject stars and artifacts. Because 
the HST point spread function is simply the diffraction limit of 
the telescope rather than depending on variable seeing, a simple 
star-galaxy separation can be achieved by plotting the SExtractor 
full width at half-maximum (FWHM) parameter versus object 
magnitude. The stars stand out as a narrow locus of objects all with 
approximately the same FWHM, up until the point where they 
saturate (at m < 18 — 18.5). For brighter stars, the locus moves 
to larger values of the FWHM while still remaining relatively 
narrow, making it easy to reject the stars. However, this star-galaxy 
separation method does not catch false detections due to stellar 
diffraction spikes and bleeding from saturated regions. Most of 
these objects can be still be flagged automatically because they 
are often highly elongated. Thus, to select real galaxies from the 
catalogs, we rejected all sources with: (1) FWHM< 0'.'13, (2) 
m < 18.5, and (3) (b/a) < 0.12, where a and b are the semimajor 
and semiminor axes, respectively. 

The COSMOS catalogs required additional flagging because 
the images obtained from the COSMOS ACS science archive con- 
tain bands several pixels wide along their left and right edges where 
cosmic rays are not properly cleaned. Simple spatial masks were 
sufficient to eliminate the spurious sources associated with the cos- 
mic rays. The COSMOS imag es were obtained at two fixed roll an- 
gles, separated by 180° (e.g.. lKoekemoer et al.ll2007l) . so that two 
sets of masking regions were required to flag the resulting catalogs. 



2.4 Definition of apertures 

We compute galaxy number counts in a set of apertures with radii 
of 45"/(\/2)\ where i = 0, 1, 2, 3. These aperture sizes are chosen 
to probe how localized any differences between the lens fields and 
control fields may be. The lens targets are centered on one of the 
ACS chips, and the i = aperture is roughly the largest that can 
fit on the chip without extending over the chip gap. For apertures 
with i > 3, the numbers of galaxies detected in the apertures start 
to drop to unacceptably small numbers on the bright end of the 
luminosity functions. In the case of the lens fields, the apertures 
are always centered on the lensing galaxy, while for the control 
fields the apertures were laid down on regular grids to maximize the 
number of independent apertures on each field while also avoiding 
the chip edges and chip gap. These grids consist of 4, 9, 16, and 
36 apertures per pointing for aperture radii of 45", 3l'.'8, 22'.'5, and 
15'.'9, respectively. Thus, for each choice of aperture size there will 
always be 20 lens apertures, wherease the number of control field 
apertures will depend on the aperture size. There will be (257 x n) 
COSMOS and (20 x n) pure parallel apertures, where n represents 
the aperture-dependent number of grid points per pointing listed 
above. 

Figure [T] shows, as an illustration, the distribution of galaxy 
number counts in the COSMOS fields inside apertures of radius 
45", with lines marking the mean number counts and the regions 
enclosing 68% and 95% of the data. The formal errors on the mean 
are very small (nearly invisible on the plot), but small number 




F814W 

Figure 1. Number counts in the COSMOS fields, using apertures of radius 
45". There are four such apertures per COSMOS ACS field. The thick solid 
line represents the mean number counts in each bin, while the negligible 
error bars on the line show the formal error on the mean. The dashed and 
dot-dashed lines enclose 68% and 90% of the data, respectively. 

statistics significantly broaden the width of the distribution at the 
bright end. Furthermore, the 2000 s exposure time for each COS- 
MOS field leads to a turnover of the number counts past m = 24. 
Thus, in the following analysis we only consider objects with 
19 < m < 24. 

Of course, having the lens-field apertures chosen to be cen- 
tered on the lens system introduces two sources of bias. The first 
is that the lens system itself contributes to the number counts 
in the aperture. The second bias is that lensing galaxies tend to 
be massive early-type galaxies and, as such, can be expected to 
be found in locally overdens e environments (e.g., [Dressier 198^; 
IZabludoff & Mulchaevll 19981) . We control the first bias by flagging 
the lensing galaxy and all lensed images by hand in the input cata- 
logs. Furthermore, we exclude any galaxies within 2"5 of the lens 
systems to avoid any strong magnification biases associated with 
these bright galaxies. The correction for the second effect is dis- 
cussed in i]4.4| There may be some additional bias in the number 
counts in the lens fields d ue to clustering associat ed with the lensed 
background object (e.g., iFassnacht et alj[2006ah , which for these 
systems is almost always a massive galaxy hosting an active nu- 
cleus. However, we expect this bias in the counts to appear only in 
the fainter magnitude bins and mostly be washed out by the large 
number of faint galaxies along these lines of sight. 



3 FREQUENTIST ANALYSIS 

In order to do an initial comparison of the three samples, we use 
Kolmogorov-Smirnov (KS) tests on three different pairs of sam- 
ples: lens vs. COSMOS, pure-parallel vs. COSMOS, and lens vs. 
pure-parallel. These tests are conducted for each combination of 
aperture size and magnitude. Figure [2] shows examples of cumu- 
lative distributions for several combinations of aperture size and 
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Table 2. KS Test Results 









logio Pks 






m 


Lens-COSMOS 


PP-COSMOS 


Lens-PP 


45.0 


19 


-0.54 


-0.00 


-0.40 




20 


-2.37 


-0.24 


-1.89 




21 


-0.84 


-0.95 


-0.63 




22 


-0.66 


-3.22 


-0.49 




23 


-0.07 


-3.78 


-0.03 




24 


-0.81 


-0.50 


-0.61 


31.8 


19 


-0.52 


-0.00 


-0.45 




20 


-1.93 


-0.12 


-1.73 




21 


-0.94 


-1.69 


-0.83 




22 


-0.11 


-8.77 


-0.08 




23 


-0.03 


-5.60 


-0.02 




24 


-0.60 


-2.04 


-0.52 


22.5 


19 


-0.09 


-0.00 


-0.07 




20 


-0.97 


-0.00 


-0.90 




21 


-1.23 


-1.23 


-1.15 




22 


-0.44 


-5.64 


-0.40 




23 


-0.03 


-1.53 


-0.02 




24 


-0.50 


-0.06 


-0.46 


15.9 


19 


-0.03 


-0.00 


-0.03 




20 


-0.23 


-0.00 


-0.22 




21 


-0.90 


-2.83 


-0.87 




22 


-0.60 


-8.18 


-0.58 




23 


-0.30 


-6.43 


-0.29 




24 


-0.13 


-5.28 


-0.13 



Results from Kolmogorov-Smirnov tests comparing pairs of samples. Val- 
ues are the logarithms of the probabilities that the given pair of samples 
could have produced the observed D values by chance if they were drawn 
from the same distribution. The "PP" designation refers to the pure parallel 
sample. Numbers in bold are those with probabiUties lower than 0.01. 



magnitude bin. Tiiese plots can be used to estimate the KS D value 
for representative pairs of samples, as well as the sample medians. 
The results of the KS tests are given in Table [51 and reveal that for 
most of the aperture-magnitude pairs the lens fields are consistent, 
at greater than the 10% confidence level, with being drawn from 
the same distribution as the control fields; only for the comparison 
to the COSMOS sample in the m = 20 bin and the 45" aperture 
is the probability that the two samples are drawn from the same 
distribution less than 1%. However, the significance of any differ- 
ences between the lens fields and the other samples is low due to the 
small number of apertures in the lens fields. With a larger sample 
of lens targets, the differences in distributions may become more 
significant. 

Somewhat surprisingly, the control samples show evidence of 
significant difference from each other, with 10 instances where the 
KS test indicates that there is less than a 1% chance that the pure- 
parallel and COSMOS samples are drawn from the same parent dis- 
tribution. Most of these low probabilities occur for the bins where 
m = 22 or 23. Furthermore, in nearly every case, it appears that 
the COSMOS fields are overdense compared to the pure parallel 
fields (e.g.. Figure |2};). While this is an unexpected result, since 
both the COSMOS and the pure-parallel fields were chosen to be 
"fair" representations of the Universe, it is perfectly possible that 
the contiguous COSMOS area is not large enough to escape be- 
ing a biased line of sight through the Universe. In fact, our re- 
sults are consistent with analyses by the COSMOS team, which 
find that, in the magnitude range 22 < i < 23, the field cho- 
sen for the COSMOS observations has higher clustering ampli- 
tudes than those found in surveys of other fiel ds jMcCracken et al.l 
l2007h . Also, the COSMOS weak lensing maps jMassev et al.l2007l ; 



iLeauthaud et al.ll2007 ) show m ore structure than seen typical sim- 
ulation fields I Faure et al]|2009l) . In this case the pure parallel sam- 
ple, albeit small, appears to provide a better indication of expected 
number counts in images of this depth. 



4 BAYESIAN ANALYSIS 

To objectively compare the number counts in the lens and control 
fields, and assess whether they are consistent, we also conduct a 
Bayesian analysis of the number count distributions. 



4.1 Poisson Fluctuations 

The first effect that needs to be considered in the Bayesian anal- 
ysis is that of counting statistics. That is, given an underlying ex- 
pectation value, /i, for the number of galaxies in a given aperture 
and magnitude bin, what is the likelihood function for the observed 
number counts, Ni, for each field i? This is simply the Poisson 
probability function 



(1) 



where P{Ni\iJ,) is already normalized. In the top row of Figure|3] 
we show representative plots that include both the distributions of 
the observed COSMOS number counts (histograms) and the corre- 
sponding Poisson distributions with the same means (solid curves). 
Figure |3jt shows that the Poisson description works well for distri- 
butions with small means, in this case a bright-magnitude bin for 
the 45" aperture. In contrast, as jj, becomes large (e.g., Figure[3j:), 
the observed distribution becomes wider than the predicted one, 
suggesting that different apertures have different underlying den- 
sity fields, i.e., the same value of /i cannot be used for all apertures 
of a given size. Clearly the analysis requires another term. 

4.2 Sample ("Cosmic") Variance 

The additional term is necessary because the presence of large- 
scale structure produces field-to-field variations, i.e. sample vari- 
ance, in the expectation value of the underlying density field. To 
model this large-scale structure term, commonly known as "cosmic 
variance" we assume that the field-to-field variations, for a fixed 
aperture, within the lens and control samp les can, in each case , be 
approximated as a Gaussian random field jBardeen et al ]|l98d) for 
H > 0, i.e. 



P{fi\lio,ao) oc exp 



2ag 



(2) 



and P{n\fiQ , (To) =0 for <= 0. The integral over the probability 
function is properly normalized to unity. The values of /lo and ao 
for a given aperture size are held fixed within each ensemble of 
fields (i.e., lens, COSMOS, or pure parallel) but can vary between 
ensembles. 

Using Bayesian theory to combine the two effects gives 



P{N,\fj.o.oo) = / P{N,\^l)P{^\^lo,oo)d^l, 



(3) 



and, combining the different fields within a given sample into a 
single data set {A*';}, with 



P(Mo,aoj{iV,}) 



P{^Jio,uo)Y{,P{N,\^lQ,Go) 



(4) 
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Figure 2. Cumulative distribution of number counts for six representative aperture-magnitude pairs. The top row shows apertures with radius 45", while the 
bottom row shows apertures with radius 15'.'9. The magnitude bin used to construct each plot is shown in the bottom right corner of the plot. The horizontal 
dashed line in each plot represents a cumulative probability of 0.5 and can thus be used to find the medians of the distributions. Blue curves represent the lens 
sample, red curves are the COSMOS sample, and green curves are the pure parallel sample. 



We assume a flat prior on fio, because it must be invariant under 
shifts, and a flat prio r on 1ok((jo), be cause it must be invariant under 
multiplication (e.g., lGregorvl[200g) . Finally, we marginalize over 
(To, which to us is a nuisance parameter, and get 

P(pq\{N,}) = j P{iiQ,cjo\{N^})dao. (5) 

To obtain the median and 68% confidence contour, we construct 
a marginalized probability distribution from P[^a\{Ni}) and find 
tiie /io values corresponding to cumulative probabilities of 0.16, 0.5 
(median), and 0.84. In Fig. [4] we illustrate the above process by an 
example for the 45" aperture and for m = 22. 

4.3 Results of Bayesian Analysis 

We calculate /.lo (i.e., the most likely average number density of 
galaxies) for the lens and control fields as a function of the aperture 
size and magnitude. The results are shown in Figure[5]and listed in 
Tabled 



Figure [5] shows significant differences between the values of 
/io obtained for the two control samples. Especially in the magni- 
tude bins corresponding to 21 < m < 23, the COSMOS galaxy 
densities are systematically higher than those seen in the pure paral- 
lel fields. This result is similar to that obtained from the KS analysis 
(Table[2]and Figure|2l( and also with the analysis by the COSMOS 
team, which finds that the amount of structure in the COSMOS 
field is at the high end of the range o f variations produced by sam- 
ple variance jMcCracken et al.ll2007l) . 

It is also instructive to plot the results in terms of differences 
with respect to one of the samples. For this exercise, the fiducial 
sample is set to the pure-parallel sample because the COSMOS 
field appears to be biased high. Figure |6] shows the resulting off- 
sets in /io. Two trends are seen in the plot: (1) For all apertures, the 
COSMOS-field values of /io are higher than the pure-parallel val- 
ues in the range 21 < m < 23, often at high significance, and (2) 
the lens fields often appear overdense compared to the pure-parallel 
fields, but the uncertainties on the lens values are so large that very 
few of the offsets are significantly different from zero. This result 
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Figure 3. Histograms of number counts in the 45"aperture for two different samples: COSMOS (a-c) and lens fields (d-f). Sliown are tlie magnitude 
bins m = 20, 22, 24 from right to left in each row. In both rows, the histograms represent the observed data while the solid curves represent a Poisson 
distribution with the same mean as the COSMOS data. Therefore, the Poisson distributions for the lens fields may not have the same means as the observed 
lens distributions. Whereas for bright galaxies the number-count variance can be explained by Poisson fluctuations, at fainter magnitudes (i.e. higher number 
counts; see panels b and c) the effect of sample variance becomes apparent. 



is consistent with our frequentist analysis of the data (§3). Only for 
the smallest aperture (15'.' 9) do we see a significant trend for the 
lens fields, with the three bins at m < 22 all being > Icr higher 
than the pure parallel values. Clearly, a much larger lens sample is 
needed. 



Overall, we can conclude that over a wide range of apertures 
(< 45"), the difference in the number of galaxies between lens and 
control field is less than ~6 galaxies. Although this can be frac- 
tionally large, it clearly shows that on average only a few galaxies 
determine the difference between lens fields and non-lens fields in 
typical observations. For the smallest aperture of 15'.'9, the differ- 
ence is typically < 1 galaxy at the 68% confidence level in any 
given magnitude bin. Since these galaxies have most influence on 
the lens model, taking these galaxies in to account has most impact 
on the lens model. We thus conclude that, given the uncertainties 
imposed by the size of our current lens samples, LSS variations do 
not make lens fields significantly overdense compared to non-lens 
fields. Only if mass does not follow light on average for the LSS 
could this conclusion be circumvented. 



4.4 Correlation-function Corrections 

In the previous subsection we saw that over the entire range of 
aperture sizes and magnitude limits, the average difference in the 
number of galaxies seen in the lens and control fields is typically 
< 6 for aperture sizes of 45", integrated over all magnitudes (e.g.. 
Fig. [6}. Strangely, even though the pointings in both sets control 
fields are random (i.e. not centered on a bright galaxy), in the case 
of the COSMOS fields one often finds more galaxies than in the 
lens fields. This might be due to the fact that, unlike either the lens 
or pure-parallel fields, the COSMOS region is a contiguous field 
that may not be a representative sample of the Universe (i.e., it may 
be overdense). 

Compared to the pure parallel fields, the lens fields are often 
slightly overdense, albeit at low significance. Based on the expecta- 
tion values in Table 3, we calculate that the integrated (over magni- 
tude) excess numbers of galaxies in the lens fields compared to the 
pure-parallel fields are 6.1±3.7, 2.2±3.8, 0.7±2.8 and 2.9±1.9, 
respectively, for the 45", 3l'.'8, 22'.'5 and 15'.'9 apertures. This ex- 
cess is not unexpected because the massive lens galaxies typically 
reside in locally overdense regions and, in principle, one would ex- 
pect more galaxies in their neighborhood. To quantify the effect 
of local overdensities, we use the two-point correlation function 
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Figure 4. An illustration of the Bayesian derivation of the most likely average number of galaxies in the lens and control fields. Shown is the result for the 
largest (45") aperture and m = 22 magnitude bin. Shown in the lower-left panel are the (normalized) number counts of lens (green) and control fields (black; 
COSMOS in this case) as function of the number of galaxies in those fields. The upper two panels show the likelihood of [fiQ, tro) for the lens (right) and 
control (left) fields, as determined from equation 4. The lower-right panel shows the marginalized probability function of for the lens (broad function) and 
control (narrow function) fields. One notices from these curves that the most likely values of for the lens and control field, in this illustration, are consistent. 



of iMcCracken et alj iimh and find predicted excesses (integrated 
over magnitude) of 2.0, 1.3, 0.9 and 0.6 galaxies, respectively, for 
the same apertures. After correcting for galaxy clustering, the dif- 
ference between the lens and pure-parallel (or COSMOS) fields be- 
comes insignificant, with differences only at the ~l-(7 level. We 
thus conclude that any line-of-sight overdensities in the lens fields, 
already being marginal, become insignificant when the local over- 
densities around the lens galaxies are accounted for. 

To test the validity of the two-point correlation correction, we 
compared two sets of number counts derived from COSMOS. The 
first set is the one that we have used as the the first control sam- 
ple, with number counts computed in grids of apertures placed 
on each COSMOS field. This set should be considered to repre- 
sent "random" lines of sight through the COSMOS fields, since 
there is nothing special about the location of the aperture centers. 
In contrast, the apertures in the second set are each centered on a 
bright (m < 20.5) galaxy found in the COSMOS area. The bright- 
galaxy sample contains 1801 apertures. In Figure|7] we plot (1) the 
COSMOS "random" number counts, (2) the bright-galaxy number 



counts, and (3) the sum of the "random" number counts and the 
two-point correlation correction. The sum of the random number 
counts with the corrections from iMcCracken et alj ilOOH ) are in 
excellent agreement with the number counts in apertures centered 
on bright galaxies. 

Hence, on average, lines of sight to lens galaxies appear to be 
overdense only because of their local overdensities and not due to 
random structures along the lines of sight. Once these local over- 
densities, and any obvious clusters along the line of sight, have been 
included in the lens models, the determination of global values of 
cosmological parameters by averaging over lens samples should be 
unbiased. 



5 COMMENTS ON INDIVIDUAL LENS SYSTEMS 

Of course, the averages can hide quite a large variation from lens 
system to lens system. Fig.[8]shows the difference between the cu- 
mulative number counts of the lens sample and the COSMOS mean 
cumulative counts for the 45" aperture. For the bright galaxies 
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Figure 5. Plots of vs. aperture size for each magnitude bin. In each plot, the blue triangles represent the lens sample, the red open squares represent the 
COSMOS sample, and the green open circles represent the pure-parallel sample. The value of /ig is determined for each aperture independently. Thus, for 
small numbers of galaxies, such as in the m = 19 bin, it is possible to have a larger fitted value of fiQ in a smaller aperture than in the next larger aperture. 
Note that the COSMOS and pure-parallel points are formally inconsistent in several of the plots (e.g., m=21, 22, and 23). 



(m < 22), the lens fields with the largest overdensities are, in order 
starting with the most overdense, SDSS J1004-I-4112, B1608-I-656, 
and B2108-I-213. If considering all galaxies with m < 24, the 
three most overdense fields are SDSS J1004-I-4112, B 1608-1-656, 
and B0218-I-357 (the three highest in Fig.[9}. All of these fields are 
outside the region enclosing 90% of the COSMOS data. The over- 
densities for SDSS J1004-I-4112 and B2108-I-213 are not particu- 
larly surprising because each of th ese systems is known to be asso- 
ciated with a cluster or rich group jOguri et al.l2004l : lMcKean et al.l 
|2009|) . However, neither B1608-H656 nor 80218-1-357 appear to be 



physically associated with such massive concentrations of galax- 
ies. Spectroscopic observations of the B 1608-1-656 field have re- 
vealed multiple group- sized associations along the line of sight 
jpassnacht etal]|2006bl) . but the B02 18-1-357 field will have to be 
examined more cl osely in future analyses. Because B0218-I-357 
iBiggs etal.lll999l) , SD SS J1004-H4112 jFghl meister et al. 20o3, 
I2OO8I) , and B 1608-1-656 jpassnacht et"aLlll999i, i2002h are systems 
for which time delays have been measured, these overdensities with 
respect to the mean line of sight must be included in any analy- 
sis to determine Ho from these systems. Fig. [8] also reveals fields 
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Figure 6. Results of Bayesian analysis, showing the difference between the values of no calculated for the pure parallel fields (green circles) and those 
obtained for the lens (blue triangles) and COSMOS fields (red squares) fields. Note the significant displacement between the pure parallel fields and the 
COSMOS fields, particularly for 21 < m < 23. 



that are underdense with respect to the mean. The most under- 
dense in bright galaxies are SDSS 1226-0006, J0816+5003, and 
B0850+054. Although time delays have not yet been measured for 
any of these systems, similar care should be taken in future Ho 
analyses that are based on the lens systems individually. 



6 INTERPRETATION & CONCLUSIONS 

To assess whether strong gravitational lenses are preferentially 
found along overdense lines of sight, we have devised a straight- 
forward Bayesian statistical number-count test which is conserva- 
tive and robust. It requires only the number counts of galaxies as 
function of magnitude inside apertures of different sizes centered 
on the lenses and on control fields. We have applied this method to 
a sample 20 lenses with F814W ACS images and control samples 
from the COSMOS and pure-parallel programs. 

Our hypothesis is that ;/ gravitational lenses are found along 
highly overdense lines of sight compared to random pointings on 
the sky, either due to structure along the line of sight or overdensi- 
ties associated with the lens galaxies themselves, then the number 
of galaxies within apertures centered on the lenses should show on 



average more galaxies than similar apertures in the control fields. 
This approach is conservative in that it does not require redshifts 
for the galaxies (hence for a given magnitude, galaxies over a wide 
redshift range contribute to the number counts), but only that the 
ratio between mass and light integrated over the redshift cone is 
close to constant. Thus, the number of galaxies can be used as a 
proxy for mass; more galaxies on average implies more mass along 
the line of sight. 

More precise statistics can be constructed if the red- 
shifts, galaxy types, etc. are known, and the galaxy masses 
are derived th rough scaling relat ions such as the Tully- 
Fisher relation I Tullv & Fished 1977) or the f undamental plane 
jPiorgovski & Davisll 19871 : Dressier et alj|l987l) . However, the re- 
sults of such models quickly become model-dependent and prone 
to systematic errors. Thus, although such approaches may show a 
more significant result than ours, a rejection of the hypothesis by 
our approach has the advantage of being a robust result that is less 
dependent on model assumptions. 

Having applied our number-count comparison to the selected 
lens and control fields we find the following results: 

(i) All distribution functions of number counts in the lens and 
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Figure 7. Comparison of number counts between the COSMOS control sample (red) and apertures centered on bright galaxies in COSMOS (green). The 
dashed line represents the sum of the control sample number counts and the two-point congelation correction factors. 



control fields are well described by a combination of a Gaussian 
random field for the sample variance (i.e., the underlying density 
field varies from field to field) and Poisson statistics in the number 
counts. This defines our underlying Bayesian statistical model. 

(ii) In the three largest apertures, with radii of 45", 38'.' 1 and 
22'.'5 (steps of 2 in area), we find no significant difference in the 
number counts between the lens and control fields (Fig. [6] and Ta- 
ble |3]l. We emphasize that this does not presuppose that there is 
no overdensity, just that our robust statistics do not require it. We 
do see a mildly significant overdensity in the ?n = 19 bin in the 
22'.' 5 aperture, which might hint that indeed there are more bright 
galaxies near the lines of sight to lens galaxies. 

(iii) The smallest aperture (15'.' 9) centered on the lenses, how- 
ever, does show significantly more galaxies in the magnitude range 
m = 19 - 21 than either the COSMOS or pure parallel fields. This 
is not unexpected since massive lens galaxies live in overdense re- 
gions (their two-point correlation is strong). 

(iv) When we correct for the effect of the two -point correlation, 
using the results from COSMOS by McCracken e t al. (2007), we 
find that a significant part of the differences between number counts 
in the lens and COSMOS fields can be accounted for by the fact that 
massive lens galaxies live in locally overdense regions, as expected. 
The remaining differences lie within the ~ l-o errors on the num- 



ber counts, even though the lens fields still show slightly higher 
number counts in the range 19<m<21in the smallest aperture. 
These differences coul d either be due to small o v erdensities along 



the line of sight (e.g.. [ Fassnacht & Lubiri 
|2006l ; iFassnacht et aljFooebl ; lAuger et al. 



2002 : iMomcheva et al.l 
20071) . or due to a too 



simplistic correction of the number counts for the COSMOS two- 
point correlation function. 

(v) On average, the excess numbers of galaxies along the lines of 
sight to the lensing galaxies, integrated over the magnitude range 
that we have explored in this paper, are small (e.g., < 6 galaxies 
for the largest aperture that we examined). These excesses, which 
are mainly due to local overdensities associated with the lensing 
galaxy, amount to only ~4-8% of the total number counts in the 
three largest apertures. These fractional overdensities can be com- 
pared with similar numbers der ived from numerical simulations. In 
particular, iHilbert et al.l 1 2007|) found, by ra y tracing through the 
Millennium Simulation l lSpringel et aTll2005[) . that lenses do along 
biased lines of sight. The contribution of the additional mass, how- 
ever, was only a few percent of the total surface mass density along 
those lines of sight, in good agreement with our results. 

(vi) While the average number counts in our sample of 20 strong 
lenses does agree well with the control samples, once the effect of 
local clustering has been taken out, individual lens systems can still 
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Table 3. Estimated value of the mean underlying number counts fiQ in the 
lens, COSMOS, and pure-parallel fields. 
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Figure 8. Difference between lens sample cumulative number counts (light 
lines) and the COSMOS mean cumulative number counts (heavy solid line) 
with apertures of radius 45". The heavy dashed and dot-dashed lines repre- 
sent the 68% and 90% ranges, respectively, of the COSMOS number counts. 
The light dashed vertical line is placed solely to identify the three most 
overdense lens fields. Starting at the top of the figure and going down, the 
line encounters, in order, the curves for SDSS J1004+4112, B1608+656, 
and B2108+213, all of which fall outside the 90% range of COSMOS field 
number counts. 
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Figure 9. Distribution of number counts in the COSMOS (light solid lines) 
and lens (heavy dashed line) fields. Plotted are the integrated number counts 
for all galaxies with m < 24 inside 45" apertures 



have significantly discrepant number counts. Thus, care should be 
taken when using an individual lens system to measure Hq . 

Based on our statistical test, we can say that: Yes, lens galaxies 
live along overdense lines of sight compared to random pointings 
on the sky, but these overdensities can be nearly fully explained 
by the fact that these massive lens galaxies are formed in locally 
overdense regions. The contribution by everything else along the 
line of sight is, in our test, not significant on average. 

Overall, we conclude that in strong-gravitational lens model- 
ing one always needs to assess the effect of the local overdense 
field, but that the effect of the uncorrelated line-of-sight over- and 
underdensities should average out in large samples (at least > 20 as 
in our case). As an example, we therefore expect the global value of 
Ho from lensing to be little affected by line-of-sight effects, but it 
could become biased if the local overdensities of the fields in which 
the lenses are embedded are not accounted for in the modeling. 
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